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Abstract 

We study various models of independent particles hopping between energy 
'traps' with a density of energy barriers p{E), on a d dimensional lattice or on 
a fully connected lattice. If p{E) decays exponentially, a true dynamical phase 
transition between a high temperature 'liquid' phase and a low temperature 'ag- 
ing' phase occurs. More generally, however, one expects that for a large class of 
p{E), 'interrupted' aging effects appear at low enough temperatures, with an er- 
godic time growing faster than exponentially. The relaxation functions exhibit a 
characteristic shoulder, which can be fitted as stretched exponentials. A simple 
way of introducing interactions between the particles leads to a modified model 
with an effective diffusion constant in energy space, which we discuss in detail. 
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I. INTRODUCTION 



Many very different glass formers exhibit surprisingly similar properties. For example, a 
common experimental feature is the 'shouldering' of the relaxation laws as the temperature is 
decreased |jl]J^ . More precisely, the relaxation of the density fluctuations evolves from a simple 
Debye exponential at high temperatures (liquid) to a two-step process at lower temperature, 
where the correlation function first decays rather quickly to a 'plateau' [jS relaxation), and then 
departs from this plateau value on a much longer time scale t(T). This relaxation time t(T) 
grows extremely fast as the temperature is decreased, in any case faster than exp^. A very 

A _ 

successful description of this divergence is the Vogel-Fulcher law: Tqt{T) ~ e^-^o, where Fq 
is a microscopic time scale [Q. However, other functional forms, such as For(T) ~ e*^^-*^, give 
reasonable fits of the data When r(T) becomes of the order of the typical experimental 

time scales (say, a day), the system is conventionally called a 'glass', which is thus inherently 
out of equilibrium. 

A remarkable and popular theory of dynamical processes in supercooled liquids is the so- 
called Mode-Coupling Theory (MCT), developped by Gotze and others Starting from a 
family of schematic equations which include a non-linear, retarded feedback of the density 
fluctuations, one can show that there exists an 'ideal glass' transition temperature Tc, below 
which the correlation function does not decay to zero ('broken ergodicity'). A notable prediction 
of the theory is the existence of the two regimes (3 and a mentioned above, and a power-law 
divergence of the 'slow' time scale t(T) as {T — Tc)^'' . However, a detailed comparison with the 
experiments shows that the transition temperature Tc, if it exists, is much higher than the 
Vogel-Fulcher temperature, leaving a whole temperature interval [To,Tc] where MCT predicts 
a partial freezing (i.e. the a regime disappears) while the experimental relaxation time is still 
finite (and behaves a la Vogel-Fulcher): the evidence for a (smeared) critical temperature is 
thus not very compelling. A way to circumvent this difficulty could be, as recently proposed in 

to work deep below Tc, where an extension of the MCT could be unambiguously tested. 

However, the physical status of the MCT is not yet very clear: the MCT equations are 
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formally identical to those describing some mean field, infinite dimensional models of spin- 
glasses where the presence of quenched disorder is assumed from the start. Further work 
should clarify the interpretation of these underlying 'spin' degrees of freedom, the importance of 
finite dimensionality effects (in particular to allow the existence of activated processes, which 
are presumably absent in infinite dimension) and the justification of introducing quenched 
disorder by hand, rather than letting it be 'self-induced' by the dynamics p. 

In view of these difficulties, it is interesting to investigate in detail alternative phenomenolog- 
ical descriptions. Among quite a few other ones [|l|,p!0|l, a particularly simple picture, advocated 
by many authors over the years, is the following: each particle is in a 'cage', i.e. a potential 
well created by its neighbours, from which it can escape through thermal activation [^[TTJ-n 



As the temperature is decreased, the probability distribution of local trapping times becomes 
very broad, and this naturally leads to a two-step relaxation, stretched exponential decays, and 
Cole-Cole susceptibility spectra. However, two important facets of these 'trap' models have not 
been addressed previously in the context of glasses: 

• At low enough temperatures, the relaxation times exceeds the experimental time scales, 
and aging effects become important |[T6|-|T8|. These aging effects have been recently discussed 
within the framework of the MCT [§|15[- Special attention will be devoted to situations where 
aging is 'interrupted' [|16|,|19| beyond a finite, but very long, time scale - which is often the case 
experimentally. 

• Since a priori all particles can move, the (random) potential well trapping any one of them 
is in fact not quenched but time dependent, further enhancing the probability of moving. In 
order to understand the glass transition, one must describe how, in a self-consistent way, all 
motion ceases. This problem is similar to the one alluded to above in the context of MCT. 

It is the aim of the present paper to discuss them in detail (a preliminary account of this 
work can be found in ^^), as well as to investigate systematically other aspects of these trap 
models, such as finite time properties. We also compare the results of this approach to the 
observed properties of glasses, and emphasize its strengths and weaknesses. 
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II. MEAN-FIELD MODEL OF NON-INTERACTING TRAPS 



A. Definition of the model 



We first consider the thermal dynamics of independent particles in a space of traps charac- 



terized by a given probability distribution p{E) for the depth E G [0,+oo[ of traps [jT6[. At 
temperature T = (3~^, each particle may escape from its trap of depth E with rate Tqc'^^ per 
unit time. In the simplest 'mean-field' version, the particle chooses a new trap of depth E' with 
probability p{E'), but with no reference to any spatial structure: this will be considered below. 
The probability density P{E,t) to be in a trap of depth E at time t therefore evolves in time 
according to the following equation 

^^^^ = -e-^^ P{E, t) + u(t)p(E) (1) 

where uj{t) = /q°° dE e~^^ P{E,t) is the average hopping rate at time t. This equation has to 
be supplemented by some initial condition P{E,t = 0) = Pq{E). 

A natural definition of the (two-time) correlation function within such a model is the fol- 
lowing: 

C{tu, + t, t^) = U{t^ + t, tu,) + qo[l - n(t^ + t, tu,)] (2) 
where n(t^ + t, t^) is the probability that the particle has not changed trap between t^ and 

tyj ~\~ t. 

n(t^ + t, t^) = £ dE P{E, t^) exp[- (roe-^^) t] (3) 

go is a certain number measuring the correlation between traps, which we shall set to zero for 
simplicity. Eq. assumes that no dynamics take place 'inside a trap' since C = 1 until the 
particle has jumped out of the trap. If a stationary distribution can be reached at long times, 
the correlation function becomes time translation invariant and reads: 

/■oo , . 

lim C{t^ + t, t^) = Ceq(t) = / dE Peq(^) exp[- (roe-'^^) t] (4) 
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B. Existence of a stationary distribution 



There exists a normalisahle stationary distribution Peq{E) at temperature T = (3 ^ only if 

1 



/ dE e+^^ piE) 
Jo 



> 



(5) 



in which case 



(6) 



In particular, Peq{E) always exists at infinite temperature (/3 = 0) and is equal to the 'bare' 
distribution p{E). The transition between the high temperature phase (T > Tq) where Peq{E) 
exists, and the low temperature phase T < Tq where Peq{E) ceases to be normalisable, takes 
place at the temperature Tq defined by 



^ = Po= lim 



E 



(7) 



We may therefore distinguish the three cases 

• Tq = if p{E) decays faster than exponentially at large E. 

• To = CX3 if p{E) decays slower than exponentially at large E. 

• To is finite in the interesting case where p{E) decays exponentially as 6"^^°^ at large E. This 



exponential form is suggested by the mean-field replica theory of spin-glasses |22|, the random 
energy model ||21[, and by phenomenological arguments in the context of glasses |jl3[. It has also 
been found exactly in the one-dimensional 'random force' model ||2^ and was used to interpret 
anomalies in the transport properties of amorphous conductors (see e.g. |2^). Its ubiquity in 
random systems is probably related to the fact that it is a stable extreme value distribution 

HI- 

For this exponential density of states, the trapping time r = Fq ^e'^^ is distributed as a 
power law for large r: 



Since this model exhibits a true finite temperature 'glass' transition, its study is most interest- 
ing. However, the case of a Gaussian density of states, where strictly speaking Tq = 0, exhibits 
several features very similar to those observed in real glasses, and will thus also be studied in 
detail. 



C. Relaxation towards the equilibrium distribution for T > Tq. 

At infinite temperature (/3 = 0), the solution of (|lD may be directly written since uj{t) = 1 

PiE, t) = p{E) + [P{E, 0) - p{E)] e-r°* (9) 

At finite temperature, we may analyse the solution of (|lD through a Laplace transform: 
introducing P(E, A) = dt[P{E, t) - P^^{E)]e-^^°\ one obtains: 

with the definition 

;(A) = / dE' e-"^ P{E',X) (11) 
Jo 



Solving for a>(A), one finds: 



^(A) 



TnA 



TOO ,p P{E,0) 
Jo "-^ l+Xef>E 



Jo "-^ l+XeHB 

In the interesting case of an exponential density of states p{E) oc e~^^^ , the denominator 
has a regular expansion up to the order A", where n is the integer part of x — 1 (x = 
followed by a singular term of order A^^^. Hence, if P{E,0) is sharply peaked (e.g P{E,0) = 
5{E—Eq)), the small A expansion of P{E, A) has a A^'~^ singularity, indicating that the difference 
P{E,t) — Peq{E) decays for large times as 

P(E,t)-Peq(^)^OC r(^-l) (13) 
t — >oo 

This decay is, not surprisingly, precisely the same as that of the correlation function Cgq defined 
by Eq. (i a 
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(12) 



C,^{t) = / dE a.eq(/5)e+^^ /3oe-^°^ exp[- (roe"^^) t] ^ ^ u;,^{P)T{x) {Tot)-^'^~'^ (14) 
JO ^ ^ 

Note that the relaxation time t(T), defined as the time after which the correlation has decayed 
down to a certain value c, diverges a la Vogel-Fulcher, i.e. 

1 



Another usual definition of the 'terminal' time, through ti = /q°° dt Ceq, would however lead to 
a power-law divergence at a higher temperature 2Tq 

r, = ^^dmrdEe^^^-^^^- = ^ (-^] (16) 

To Jo T>2To To \T-2Tq) ^ ' 

The existence of two characteristic temperatures in this model was recently emphasized in 
the context of glasses by Odagaki [Jl3[. If p{E) decays faster than exponentially (say as 
exp(— ii^^'^'')), the relaxation towards equilibrium is faster than any power-law, since the small 
A expansion of P{E, A) is regular. However, slower than exponential relaxation (such as 
exp — (logt)^'*"'^) is expected. More generally, the long time behavioiu' of P{E,t) — Pcq{E) 
is the same as that of Ceq for large t. 

D. Non normalisability and aging in the low temperature phase. 

When there is no equilibrium distribution (|^), one expects that the dynamics will never 
become stationary [|1^,^,^,Q , and this leads to aging effects. One may look at large time for 
a scaling solution for the equation (|1|). We introduce the natural dimensionless scaling variable 
u = s-i^d a dimensionless function normalized as /q°° du (j){u) = 1 and try to find the 
asymptotic distribution function through the form 

P{E,t) ^ (3u(j){u) (17) 



The resulting equation for is 



The left-handside presents a t independent limit as t — > oo only if the integral on the left- 
handside diverges in such a way that it compensates exactly the decay of p |^l2iil2M^ as t — > cxd. 

Consider for example the case of an exponential density of states p{E) = (3q e"^"^, for which 
there is no equilibrium distribution when x = ^ e]0,1[. 

In order to get some sensible limit of equation ([T8|) as t ^ oo, the function must present 
the singularity 

<P{u) ^ 7 (19) 

n— >0 

where a is some normalisation constant. The solution of the resulting equation 

+ (m - 1)0(m) = -7 (20) 

(Jj (Jj 

satisfying the normalization condition du (l){u) = 1 reads 

_1 1 

0(^) = ^-e ufudyy^-'ey (21) 
vri ix) u Jo 



The consequence of the scaling form (1?) on the correlation function C or 11 defined as the 
probability for the particle to remain in the same trap during the interval + 1] (see Eq. 

(D) 

n(t^ + t, t^) = r dE P{E, t^) e-(^'>'~'''y (22) 
Jo 

is the aging behavior 

1 / t 



n(t^, t^ + t)- r du (f){u) e u\tu,J (23) 
Jo 



where the two times only appear now through the combination (^j- This must be contrasted 
with the situation prevailing for T > Tq, where as soon as ^ t(T), P{E, tw) ceases to depend 
on tw and the correlation function only depends on the time difference t, and is plotted in Fig. 
1-a for different values of ^ (We have included a 'fast' relaxation inside each trap - seen as a 
simple harmonic well, in order to mimic the j3 regime of glasses). 
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For the case of the exponential density of states, using equation (pTf) we recover the exphcit 
formula derived in through another approach 



j_] av [I — vj 



n(t^ + t, t^) ~ — — / (^) dv (1 - vf-' v-^ (24) 
with the asymptotic behaviors 



l-x 

-] 



U{tu, + t, t^) ~ 1 - sinc[7r(l -x)][^) for t < (25) 



n(t^ + 1, ty^) ~ sinc[7rx] (—) for t tm (26) 

where we have introduced the notation sincfwl = 

L J u 

E. Case of Gaussian density of states at low temperature : 

stretched exponentials and interrupted aging 

Let us consider the interesting case of a Gaussian density of states 
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- 



p(E) = -|-e" ^EJ (27) 

which corresponds to a situation where the energy barriers results from a (weak) interaction 
with many neighbours. Strictly speaking, there is a stationary distribution for any positive 
temperature (H) 

/E\^ {E-E*Y 
P^^{E)=u,^{f3)^-e+^^ e~\E'J ^Ue where E* = ^{(3E,) 

(28) 

However, at very low temperature T <^ Ec, there exists an approximate aging behavior on the 
time interval satisfying 

< In(rot) < ifSEcf (29) 
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but this aging phenomenon progressively disappears as time becomes large. The physical 
reason is the following: as can be seen directly from p{E), the distribution of relaxation times 
r = Tq^c^^ can be written (for large t) as 

-I In(roT) 

^(r) oc r (30) 

Hence, the parameter x = ^ defined in the previous paragraphs is replaced by a slowly varying 
function x(t) = 2 ^°|^°2^ ■ (The factor 2 comes from a more careful analysis - see below). The 
above inequality thus corresponds to the case where x{t) <^ 1, where aging effects are indeed 
expected. 

More explicitly, one may look for an approximate scaling solution of (|I|) analogous to (|lj 



P{E,t)c^PuMu) (31) 

/3-E 

where u = and where the function 0^ is supposed to vary slowly in time. More precisely, 
in the resulting equation for 0t 

^■^i 2^0t \ 1 / log (rotw) \ dv ( V 



we assume that the first term containing the time derivative is negligible in comparison to 
others. We now therefore consider the simplified equation 

u^ + {u- l)Uu) - P j / (f^ j (33) 

and we will discuss the validity of the hypothesis made on the exphcit solution obtained for 

m- 

For TqI ^ 1, the function (pt must present the following asymptotic behavior 

0t(M) ~ A(t) e V/^^J (34) 

and (|33D becomes 

+ {u- l)<j)t{u) ~ A{t) e Kf^^J (35) 

(jj Uj 
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The solution reads 



1 



A{t)- e u 
u 



1 1 

u dy y^W-^ 



'logy' 







(36) 



where A{t) is determined through the normalisation condition /q°° du(j)t = 1 

]_ ^oo r+oo pzx(t) 



Ait) 



dy e 



+ 00 

dz 



+ y 



The asymptotic behaviors of (t)t{u) are respectively given by 
for u ^ oo 



for u 



(37) 

0, and by the following 



xit) 



M-i-^W e 



'logw' 



(38) 



We may now estimate the order of the term dropped in (|3 



u t 



dx{t) d(f)t 



dt dx{t) 
The comparison with other terms of ( 



u 



{(3E,y dxit) 



2 logu 



(39) 



shows that the approximation ( |33D holds at low 



temperature (PEc)"^ ^ 1, and not too large u, i.e.: 



log u 



< 1, 



Let us now explain how our analysis is compatible with the relaxation towards the equi- 



librium distribution (|28|) . At low temperature {(3Ec ^ 1), this equilibrium distribution is 
concentrated in a region of scale Ec around the central value E* = which is very large 
compared to E^. This is why during a very long time interval, the probability distribution 
P{E, t) is pushed towards larger and larger energy as if there were no equilibrium distribution. 
More precisely, as soon as T^t ^ 1, P{E,t) enters an aging-like scaling regime in the variable 



u 



Tot 



as in the case of the exponential distribution below the transition (T < To). From the 



above remark concerning the shape of ^'(t) (Eq. |3^), one expects that the aging behaviour 
disappears as x{t) gets close to the value 1. Then the normalisability of Pcq{E) becomes notice- 
able, the probability distribution P{E, t) begins to stabilize, and the scaling variable u = 
looses its meaning. In the limit x{t) — > oo, the function 0j becomes concentrated in the region 
M — >• 0, and the asymptotic behavior (|3^ precisely corresponds to the equilibrium distribution 
(|28|) when one replaces u by 
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We now turn to the correlation function 11 defined in (|22D . For TqI^ ^ 1, one has in terms 
of the function (pi^ (u) of equation ( ^6l) 

U{t^ + t,t^) - r du <PtM e « Vt^y (40) 
Jo 

Here, contrary to the real aging behavior (p3[) , n(t^ + t, 1^) is not a function of only, since 
the function (pt^ slowly varies with the waiting time t^. In particular, using the asymptotic 
behaviors (0) and (|38D, one has respectively at short time t <^tw 



f log u ' 

n(t^ + t,t^) ~ l-A{t^) r du e 
and at large time t 3> 

A{tyj 



I 



1 _ e u\t 



(41) 



logy' 



(42) 



To get more explicit behaviors, we specialize to the interesting regime < xit^ < 1 where 
the aging-like behavior occurs. In this regime, the normalization constant ( P7D may be expanded 
according to 



1 



)roo r+oo 
Jo Jo 



z + y 



(43) 



'J d_ 



sin 7rx(t^) 



(44) 



The explicit expression of the first correction in the small parameter -jj^^ shows that the 
approximation 

sin 7ra;(t^) 



7lT{x{tyj)) 



(45) 



is vahd for -C a;(t^) ^ 1 — 



The two asymptotic expressions (|4ll - ^2D may be analysed in the same way. At short times 
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1 d V 



1 "^^y^Xli^ \ 



(46) 



can be approximated by 



/ t \ 

n{t^ + t,t^) 1 - sinc[7r(l - — (47) 



for -C x{tu,) ^ 1 — and for ^ not exponentially small in {PE^y. 
At large times, 



1 d V 



n(t. + M.)^^ iMe"U^c^-(^. 



T{x{t^)){-' 



(48) 



may be reduced to 



/ t \ ^^(*«') 

Il{t^ + t,t^)^ sinc[7rx(t^)] — (49) 

for -C x(tw) ^ 1 — and not exponentially large in {PEcY. 

As anticipated, these expressions are very similar to those obtained for the exponential 
density of states, provided one defines an effective time dependent parameter x{t). [Actually, 
when I In^l ^ 1, a slighly more accurate expression of (^^91) is obtained by replacing the 
exponent x{tw) by in agreement with the naive interpretation of Eq.(pOD] 

Let us now concentrate on times such that xit^) > 1: aging effects cease, and the stationary 
dynamics resume. The resulting correlation function 

Ce^{t) = l rfE a;eq(/3) e^^-^e-(^) exp[- (loe-'^^) t] (50) 

is plotted in Fig. 1-b and displays very interesting features. To a very good approximation, 
one may replace the effect of the last exponential by a simple cut-off: 

f3E f_ 

One thus finds the following behaviour: 



Ceq(t) ~ I dE W,q(/3) e"^ ^-6-^^ ) (51) 



aCeq(t) 1 \ PE, 



d\n{Tot) ^PE, 
13 



(52) 



where t(T) = exp is the time scale associated with the a regime. As mentioned in 

the introduction, the experimental data on r(T) can indeed be fitted by such a law [^. Note 
however that the 'terminal' time scale, defined as ti = /q°° C^qdt, and extracted, in supercooled 
liquids, from viscosity measurements, is given by ti ~ jr-el(/3^c)^ temperature [|. 

Eq. (|52|) shows that the correlation function measured for different temperatures approxi- 
mately fall onto a master curve if one plots them as a function of provided the variation 



of l3Ec (compared to that of t(T)) can be neglected. Furthermore, as emphasized in the 



shape of the relaxation function given by Eq. (|52D is actually very close, in the region where t 
and t(T) are not too different, to a stretched exponential^ Ceq{t) = exp — (7^^) , provided 
one identifies a{T) ~ (1 + I^)"! |2^. 

These features are again reminiscent of the experimental ones, where an approximate scaling 
of Ccq{t) for various temperatures can indeed be achieved, and where a stretched exponential 
(or Kohlrausch) form for the decay in the a regime is quite often proposed (see the in [Q). 
Interestingly, the exponent a extracted from experimental fits has a tendency to decrease as 
the temperature is decreased (see, e.g.. Fig. 16 of [^). 

Finally, let us stress that many of the results of the present section are not restricted to 
the Gaussian form of the density of states, and would still hold if p{E) oc exp — (-^)^+'^, u > 
with {(5EcY is replaced by {(3EcY^'^' (up to v dependent prefactors). Such a generalized model 
has been recently discussed in the context of glasses in p9[. An interesting point to notice 
is that the variance of the log-normal form of ^^"-^^ is related to its maximum through 



^ Note however that the long time behaviour of Ceq(t) is well approximated by Ceq(t) a (Jj-)^ with 
H = (-^)^ log Tot and T2 = Tq ^ exp(/3i?c)^) showing that there is actually quite a bit of freedom to 
choose the rescaling time factor t{T) without affecting too much the quality of the scaling. 

^Note: the stretched exponential exponent is usually called j3, although it describes the a relaxation 
in glasses, and should furthermore not be confused with the inverse temperature. We have thus 
proposed in to call it a. 
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cr^ = Mnr(T). Such a relation was also obtained by Souletie using different arguments |l30[| , 
and is not incompatible with experimental data. 

The conclusion of this section is that stretched exponential decay and scaling of the relax- 
ation curves can occur in the absence of any criticality or cooperative effects, and arises even 
when the distribution of local energy barriers decays quite fast for large E. An interesting 
consequence is the appearance of 'interrupted aging' effects for low temperatures, with a cor- 
relation function described by Eq. (|T7|Ji9|) . Note that this correlation nearly scales as but 



with a systematic bias similar to the one discussed in |19|. 
The above model is however oversimplified since: 

- no spatial structure is included and 

- the particles are independent. 

The next two sections aim at discussing, at least partially, these aspects. 

III. MODEL OF TRAPS IN EUCLIDEAN SPACE 

A. Continuous-time random walks 

We now consider that the traps live on a d-dimensional hypercubic lattice. The particle 
performs a usual random walk on this lattice, but has to wait for a certain trapping-time r 
before each jump [^]. For a trap of depth E, the distribution of the trapping-time r reads 



vl/^(r) = To e-^^ e~('^'^ (53) 

If p{E) denotes the probability distribution of the depth E of the traps, the particle performs 
a random walk among traps with a distribution of trapping times 



oo 



^(r) = / dE p{E) ^E{r) (54) 
Jo 

Note that we consider here a model with annealed disorder, where the energy E for a given site 
changes at each visit of the particle. This is justified for dimensions d > 2, where the particle 
rarely visits twice the same site. The model that we consider is thus the well known 'continuous 

15 



time random walk' (CTRW) first introduced by MontroU and Scher, and discussed many times 



31| . However, the non stationary properties of tlie low temperature phase (when it exists) has, 
to our knowledge, never been considered - only the case where the waiting time is zero has 
been investigated. We shall thus focus on the following natural correlation function: 

t + t^, t^) =< e^9^-(-l*+*-)--^(*»)) > (55) 
in particular in the case where there is no equilibrium distribution for (|ID. 

B. Correlation function C{q,t,tw = 0) 

We shall first recall the result for a vanishing waiting time = 0, which is the case usually 
considered |31|. One has by definition: 



C(g, t, 0) =< e'^-("^(*)-"^(°)) >= J2 e''" ^(^' ^|0' 0) (56) 

r 

where the summation is over all sites of the lattice. The probability density may be decomposed 
according to the number of jumps as 

OD 

P(f,t|0,0) = QNi^RtiN) (57) 

N=0 

where Qn{''^ denotes the probability to be at site r after jumps for a usual random walk on 
the lattice. Rt{N) denotes the probability to perform jumps in time t and reads in terms of 
the distribution \E'(t) of trapping time and Heaviside function 6 

Rt{N) = ( n dT,^{T,)\ 0\t-J2J^je\^2T,-tj (58) 

or more explicitly in terms of the probability density \l'Ar(r) of the sum r = J2f=iTi oi N 
independent trapping times 

Rt{N) = f^dT(^N{T) - ^N+i{r)j (59) 

The correlation function may then also be decomposed according to the number of jumps as 
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C(g,t,0)= J2MN)QN{q) (60) 

N=0 

where Qn{q) denotes the Fourier transform of a usual random walk on a d-dimensional hyper- 
cubic lattice of lattice spacing a starting from site 

/I .AT 

QN{q) = E e^'"^*^ ^iv(r-) = h E (61) 
It is convenient to introduce the Laplace transform 

(7(5, A, 0) = / dt e-^'C{q, 0) = E QNiq)RxiN) (62) 
•^0 JV=0 

where 

roo 1 / \ iV r 



~ 1 / ~ 



1 -^(A) 



(63) 



in terms of the Laplace transform ^'(A) of the distribution ^(r) of trapping time 



oo 



*(A) = / dr e-^^ *(t) (64) 







We obtain finally 



1 1 - '^(X) 

^i^^ 0) = A TT^ ^ ^^^^ 

l-*(^)(^Ecos(g^a)j 



C. Correlation function C{q,t + tw,tw) 

We now turn to the case of arbitrary waiting time where the correlation function reads 
now: 

POO 

C{q, t + t^, t^) = E E e*'^^-^^"'*'"^ J dE^ r{f, t + t^\f^, E^, t^)V{f^, E^, t^\0, 0) (66) 
It is convenient to reorganize this expression as follows 

C{q,t + t^,t^) = J dE^Y. V{f^,E^,t^\Q,Q)^ e'i<r-r^) r{f,t + t^\f^,E^,t^) (67) 
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and to decompose the conditional probability V{r, t + t^|r\„, Eyj, t^) into 

— fFne"^^™^ t 1 
V{f,t + Uf^,E^,t^) =6r-r^ e V ; +y ^p(f,t_^|f^ + ^,0) (68) 

The first term takes into account the probability to stay in the trap of energy during 
the time interval + 1]. The second term describes the probability to jump at the time 

tw + t out of the trap of energy i?^ to go to one of the 2d neighbours fv, + e of f\„, and where 
the 'aging' process starts anew. One has therefore 

E ^^''''^ t + Ur^, E^, t^) = e ^ (69) 



r 



Noticing that 



dr-^EM E cos(g^a))c(g, t - r, = 0) (70) 



E V{r^, E^, t^\0, 0) = P{E^, t^) (71) 



is simply the probability distribution studied in the previous section (|T]), and that 







dE^P{E^,t^)e V 0^ )=U{t + t^,t^) (72) 



is the correlation function (^) introduced previously, we obtain finally 

C{q, t + t^, t^) = U{t + t^, t^) - H E cos(g^a) j J dr C{q, t - r, = 0) ^n(r + t^, t^) 

(73) 

This convolution product leads to introduce the Laplace transforms 

C{q, A, t^) = / dt e~^' C{q, t + t^, t^) (74) 

and 

n(A, t^) = I dt e-^* U{t + t^, t^) (75) 

to get the simpler relation 
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Equation (^5]) gives finally 



1 - ^ ( ^ E cos{q,a))c{q, A, 0)1 + f ^ E cos(g^a) )C{q, A, 0) (76) 



C{q,X,t^) =fl{X,U 



1 - I - ^ cos(g^a) 



1 -^(A) 



^1 E cos(g^a) 



(77) 



^(A) 



Another quantity of interest in glasses is the susceptibility, which is related to the above 
correlation function through: 



x(g, 00, t^) = 1 + iuj / dt e*"^ C{q, t^, + t,tyj) = 1 + iu C{q, -iuj, t^,) 



(78) 



(see however M for a discussion of the Fluctuation-Dissipation theorem in this context). 



D. Case of exponential density of states at low temperature 



For the exponential density of states p{E) = /3o exp(— /^o-E), there is no equilibrium distri- 
bution for the process (|I]) when x = < 1. The correlation function 11 is known in this case 
(|2^) and its Laplace transform reads 



n(A,t^) ^ — / du e 
Jo 



IT 



with the asymptotic behaviors 



fl{X,t^) 



1 (At, 



At„,«i r(i + x) A 



dv (1 — v)^ V 



(79) 



(80) 



n(A,t. 



1 

At„,>i A 



1 1 



r(x) (At^)i- 



The distribution of trapping times 



HdE p{E) Fo e-^^ e-(^° 
Jo 



xFfi 



(ror 



For 



dy y^'e 



(81) 



(82) 



presents the slow algebraic decay 
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^ , , xT(l + x) , , 

^ r ~ ^ , ' (83) 

and its Laplace transform presents therefore the non-analytic behavior 

/■oo 1 / A \ 
^(A) = / dr e-^^ ^(r) ~ 1 — (84) 



Moreover, one has for wavevectors such that ^ is much larger than the lattice spacing a 

1 



9 

d 



Y^oos{q,a)) ^ I- -/ + ... (85) 
It is convenient to introduce the time scale tq 

(r„t,r = ^ (86) 

corresponding to the typical time needed by the particle to spread over a region of size 1/q. 
There are therefore three time scales in ( [77D : A^^, tg, t^. We are interested in the region where 
all three are much bigger than the microscopic time scale Tq ^ : A~^ ^ (fT^-p^) ,tq':^^ (p5| ) 
and t«, ^ (fT^D- Still we have to distinguish various regimes in (jn\). We give the asymptotic 



results for the correlation function C{q,t + t^,t^) after the inversion of the Laplace transform 



t 



l-x 



C{q, t + t^, t^) ~ 1 - sinc[7r(l -x)][ — ] for < t < (87) 



C{q,t + tw, tw) — sinc[7rx] i—\ for tq <^ <^ t (88) 

sine TTx t 

C{q,t + t^,t^) fort<t^ andt^t, 

T[x) t^tir^ 



(89) 



, _ sincfTTxl ( 
C{q,t + t^,t^) - 1 - M - fort^<t< 

r(l + a;) \tq) 



tq (90) 



C(g,t + t^,t^) ~r(l + a;) I - ) for < < t (91) 



^1, 
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There are three interesting points to notice: 

• For qa very large, such that tq <C t,t^, we find the same asymptotic behaviors as for 



n(t + tyj) (^-26) in the previous section. Physically, this means that as soon as the particle 
has jumped once, the rapidly oscillating correlation function averages to zero. Hence, only the 
particles which have not yet moved contribute to the correlation. 

• There are two regimes where the correlation function behaves similarly to a stretched 
exponential at small times, when tq <^ t <^ t.^ or tyj <^ t <^ tq. The exponent a of this 
stretched exponential is however different in both cases: it is equal to a = a; when <^t <^tq, 
and equal to a = 1 — a; in the other case. 

One should however stress once more that the present trap model is unable to explain why 
there should be a relation between the exponent a describing the decay of the correlation in 
the a regime which we discuss here, and the shape of the relaxation in the short-time {(3) 
regime, which corresponds to intra-trap dynamics. That such a link exists is one of the major 
prediction of the MOT, which should also exist deep in the glass phase, where aging effects 
similar to those discussed here are present. 

• The frequency dependent susceptibility defined by Eq. ([78|) behaves in a Cole-Cole fash- 
ion: 



for uj <^ Fq, qa <^ 1. This has been emphasized in e.g. [0. However, in the aging regime 
ojtyj ^ 1, the behaviour of x is given by: 



A similar expression was obtained in the context of spin-glasses in pSj . 



IV. MODEL OF INTERACTING PARTICLES IN TRAPS 
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A. Definition of the model 



In the above sections, we have considered models for which the motion of a given particle 
out of its trap does not affect the potential seen by the others. In order to model this effect, 
we have proposed in to add to model (|l|) a diffusion term in energy space proportional to 
the mean hopping rate itself u{t) = /q°° dE e~^^ P{E,t). The equation for P{E,t) then reads: 



1 (9P 8 

-^^-e-^^P + .mE)+.it)D- 



p{E) 



dP 



P 



dp 



(94) 



BE dE_ 

This diffusion term expresses the fact that every 'hop' induces a small change in all the 
neighbouring e's. Assuming that the transition rate is proportional to the final density of 
states, the contribution of such an effect to the master equation a priori reads : uj{t) J de' T{\e — 
e'|){P(e', t)p(e) — P{E,t)p{e')}. In the limit where the width of T(|e — e'|) is small, justified in 
a mean-field limit where the number of neighbours is large, this term reduces to the diffusion 
like term in (0), with an effective diffusion constant D proportional to the width of T. More 
general forms for this diffusion term will be discussed below. 

As before, equation ( plD has to be supplemented by some initial condition P{E,t = 0) = 
Pq{E) and by a 'hard wall' boundary condition aX E = Q 







E=0 



(95) 



to ensure the conservation of probability. 



B. Existence of a stationary distribution 



There exists a stationary distribution Peq{E) at temperature T = (3 ^ only if the equation 



MEl 
dE^ 



(96) 



piE) u,,p{E)_ 

supplemented by the boundary condition (R^) admits a normalizable solution. The discussion 
again depends on Tq = I3q^ introduced in (|^. 

• For T > Tq, the asymptotic behavior at high energy of the solution of (|96|) reads 
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(97) 



which is normahzable since j3 < j3o. 

• For T <Tq, equation ( p6D can be approximated at high energy by 



d^Pe^iE) p"{E) 



Pe^iE) 



1 
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dE"^ p{E) 

The general solution reads using arbitrary constants A and B 
Ag,iE)+Bg2iE) + ^ 



(98) 



(99) 



f"^ du 
p^{u) 



(100) 



9i{E) I dv g2{v) + g2{E) I ^ dv gi{v) 

in terms of the two independent solutions 

gi{E) = piE) and g^iE) = p{E) T 

Jo 

of the homogeneous equation. Taking into account that p{E) decays at least exponentially at 
large E for /3o < +oo, one may show that there is no normalisable solution. 

In summary, the conditions for the existence of a stationary distribution in term of Tq are 
therefore exactly the same as in the case D = (^. 



The case of an exponential density of states was considered in [0]. For completeness, we 
give the explicit form of the equilibrium distribution: 



Peq(^) = AT 



Kjz) 



(101) 



Ki,_i{zo) 

where u = ^^fr, z = zq exp and zq = ^^S. and K^, are the Bessel functions of order 
u, and lCy{x) = ^K^{u), T^{z) = ^I^{u). J\f and F are fixed by the normalisation of 
Peq{E) and the boundary condition (pSj), which lead to the following equation: 



D 



Iu-l{zo) 



1{Zq) Jzo U 



(102) 



C. Aging in the low temperature phase 



When there is no equilibrium distribution (^, we may proceed as in the case D = (section 
III D|) and look for a scaling solution of the form ([T7|) for the equation (p4|). The resulting 
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equation for the dimensionless function again only admits some non trivial limit as t — > cxo 
if the left-handside of (^) does. For example, in the case of an exponential density of states 
p{E) = (3q e~^°^, in which there is no equilibrium distribution when x = ^ g]0, 1[, the function 



must present the singularity (|T9D 



(f)(u) ~ 7 u-'-" (103) 

u— >0 

The equation for now generalizes equation (pO|) 

^^'tI + [3^"' + ^^'^1^ + [^(1 - + iu- = -7 (104) 

where A = 'yDf3^. The normalisation constant 7 is determined by the normalization 
/q°° du (j){u) = 1. Contrary to the case D = (PTf), the solution cannot be explicitly written. 
However, its asymptotic behavior at large u is easily obtained from (104) 



m) ~ ^ M-i-^ (105) 



1i— >oo 



X 



The aging behaviour (^) of the correlation function (^2]) still holds as a consequence of 
the scaling form (0). The asymptotic behaviors (|103[ - |105D of the scaling function induce the 



asymptotic expressions generalizing ( p5| - 126|) for the correlation function 



n(t^,t^ + t) ~ 1-7 7^3^ (t^) fort<t^ (106) 

n(t,„,t,„ + t) ^7 ^ (t^) fort>t^ (107) 

The presence of the diffusion term in (Q) therefore only affects the normalisation constant 7, 
but does not change the asymptotic time dependence of the correlation function Hit^.t^ + t) 
of the model (|1|). This suggests that the difference between the model considered in section 
(|11 D| ) (corresponding to = 0) and the 'annealed' model considered here is, to some extent, 
irrelevant. 
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D. Possible generalizations of model (^) 



The model (W) may in fact be generalized through the introduction of an energy dependent 



diffusion constant, i.e.: 

(108) 



^dE dE 

with the same boundary condition (pSf ) as before to ensure the conservation of probability. p{E) 
still corresponds to the equilibrium distribution at infinite temperature (/? = 0). The choice 
D{E) = for example, would correspond to an 'entropy' biased diffusion in energy space, 
with a driving force proportional to '^|^°J^ , and an effective diffusion constant independent of E. 

A careful study of the existence of a stationary solution to this equation shows that there 
are two cases. 

• If lim — log D(E) < 0, the transition takes place as before at the inverse temperature (3o 

E^oo E 

• If lim — log D{E) > 0, there exists a stationary distribution at any finite temperature. In 

-B— >c« E 

the case where D(E) = and p{E) exponential, for example, one finds that the large E 
behaviour of Peq{E) is of the form: 

PeJE) ~ Ae-'^o^ + 5e~2^°^ (109) 

showing that the correlation function now decays as Ceq(t) oc t ^ ' for large times. Note 
that the terminal time scale ti = /q°° dt Ccq(t) diverges when T < Tq, although no asymptotic 
aging effects appear in this temperature regime. 



V. CONCLUSION 



We have studied in this paper a model of particle hopping between energy 'traps' with 
an arbitrary density of energy barriers p{E). As emphasized in |1T6|J25| , P^ , the case where 
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p{E) decays exponentially is special because it leads to a true dynamical phase transition 
between a high temperature phase and a low temperature aging phase. More generally, however, 
on expects that for a large class of p{E), 'interrupted' aging effects appear at low enough 
temperatures, with an ergodic time growing faster than exponentially. It would be interesting to 



look systematically for aging effects experimentally [17]. Furthermore, the relaxation functions 
have a strongly stretched shape (see e.g. Fig. 1-b), which can be fitted as stretched exponentials. 
The case where the traps are organized on a c? dimensional lattice is slightly more involved, since 
a new time scale tq appears, which is a wave- vector dependent relaxation time. A schematic way 
of modelling the interactions between the particles, reflecting the fact that when one particle 
moves, the potential energy seen by its neighbours changes, was investigated in section 4. The 
conclusion regarding the existence of a dynamical transition was found to depend on the shape 
of the effective diffusion constant D{E). In the case where the transition survives, the role of 
this 'interaction' is irrelevant. 

In conclusion, many of the observed features of glassy systems can be accounted for within 
a picture of independent particles trapped in random potential wells, without any obvious 
collective effects [^p!T|-[T^. How deep is the link between this picture and the Mode-Coupling 
Theory is not yet clear. However, the MCT relies on the existence of a true critical point 
, with two important consequences which are beyond the grasp of 'trap' models: the strong 
link between the a and (3 regime and the critical behaviour of the plateau value in Ceq(t) as 



a/Tc — T [|T|. These predictions are however hard to test since the critical point is supposed 
to be blurred by 'activated processes'. As suggested in [Q, a crucial test of MCT could be 
performed by working deep below Tc, where 'activated processes' are frozen but where strong 
aging effects should appear. 

Appendix A 

Relaxation towards the equilibrium distribution in model ( |94D 
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In the case D > 0, equation (0) is nonlinear and cannot be solved exactly. However, we 
may introduce some approximations to study the relaxation at large time. 

Approximate description of the relaxation 

To study the relaxation towards the equilibrium distribution Peq{E), we may set 



P{E,t) = Peq{E)+T]{E,t) 



(110) 



into (|9^) , linearize the corresponding equation in the perturbation rj, and replace Pcq{E) by the 
expression Peq{E) ~ LJcq{l3)e~^^^ p{E) which happens to be a very good approximation except 
in the vicinity of E = 0. We obtain finally 
1 dr] 



Iq ot Jo 



+o;eq(/3)D 



d^E 



VP"{E) 



(111) 



We may look for the solution through its decomposition onto a relaxation spectrum 



with the condition 



v{E,t)= / dXfx{E)e~''^'' 
Jo 



POO 

/ dEfxiE) = 
Jo 



'112) 



(113) 



to insure the normalization of the probability density P{E,t) ( |110| ). The equation for fx ( |1 1 1| ) 



fx{E) = Cx 



(114) 



where 



(115) 



has to be supplemented by the boundary condition (I 



fj^^dfx dp 
P^^'^dE-f'dE 



E=0 



(116) 
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In fact, the self-consistency condition (|115|) is automatically satisfied once f\ is solution of 
(|TT1 - |TTWT6|) . 



Equations ([TT4| - |Tl3HTl6D imply that 



(A -AO r^^h(E)fy{E) = (117) 







The functions fx{E) including Peq{E) = fo{E) are therefore orthogonal to one another with 
respect to the measure . One may choose to orthonormalize this set of functions according 
to the scalar product 

dE 

ME)h,iE) = 6{\-\') (118) 



/o piE) 

The development of P{E, t) onto this basis then reads 

P{E, t) = Peq(^) + / dX a{X) h{E) e-^r°* (119) 
Jo 

where the coefficients a(A) are simply obtained through scalar product with the initial condition 
P{E,t = 0) 

«(A)= I ;^/A(i?)P(i?,0) (120) 
Nature of the relation spectrum 



To determine the relaxation spectrum, it is convenient to transform (|118|) into a usual scalar 
product 

roo rlTT' POO 

S{X - A') = X -7^ fx{E)fy{E) = du ^x{u)^Au) (121) 



through the change of variable 



dE' 

E u{E) = / ^= (122) 



or conversely 
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Eiu) 



and a change of functions 



" dv 



h{v) 



where 



h{u) 



p[E{u)] 



(123) 



ME) <ilx{u) 



fxjEju)] 
{p[E{u)])'^ 



(124) 



The equation for the new function ^'a('") ( |114D 



+ 



V{u) + 



.-f3E{u) _ 



(125) 



where 



Viu) 



15 fh'iu)\ 3h"{u) 



4 V Hu) J 2 h{u) 
has to be supplemented by the boundary condition at uq = u{E = 0) 



(126) 



d-^x 3h'(u) ^ 

— - H —'^x 

du 2 h{u) 







(127) 



and by the condition ( |113| ) 



ii 



(128) 



'-"0 [^(^)]- 

To discuss the nature of spectrum, the only important property of the potential V{u) is that 



V{u) 







(129) 



This comes from the fact that p{E) presents a rapid enough decay as — > oo to be integrable. 

At infinite temperature (/5 = 0), the constant Cx ( |115| ) on the left-hand side of ( |125|) vanishes 
( |113| ). The asymptotic form of equation ( |125| ) as m ^ oo 



d^^x ( \-l 



^xiu) 



(130) 



admits two independent oscillatory solutions for A > 1, and the boundary condition at Uq 
determines the suitable linear combination to be takken for "^x up to a normalization constant. 
However, for A < 1, only the exponentially damped solution is acceptable at infinity and the 
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boundary condition at Uq cannot be satisfied. So we finally obtain that the relaxation spectrum 
at infinite temperature [jS = 0) consists in a continuum above the gap Tq 

/oo 
dX a{X)ME) e-"r«* (131) 

to be compared with the previous result for the limiting case D = 

At finite temperature {(3 > 0), the asymptotic form of the homogeneous equation corre- 
sponding to ( |125|) in the limit u ^ oo 

admits always oscillatory solutions for A > 0. The relaxation spectrum at finite temperature 
{(3 > 0) consists therefore in a continuum starting from A = (|119D. This will produce an 



algebraic relaxation determined by the behavior of the function (px in the limit A ^ 0. 

Example of the exponential density of states 

We may apply the previous general theory to the particular case of the exponential distribu- 
tion for p{E), for which there exists a stationary distribution in the domain of high-temperature 
X = ^ > 1. The change of variables ( |122| ) 

E — > u{E) = uq where uq = — '^—^ (133) 

and the function h{u) (|123|) 



2 



h{u) = (134) 

are simple enough to give an algebraic form to the potential involved in the inhomogeneous 
Schrodinger equation ( |125| ) satisfied by the new function \E'a(m) 

d^^x 



+ 



15 1 1 f Un\x A 



^A(«) = 7T^f^l (135) 



At infinite temperature (a;eq(/5) = 1 ; x = oo), the function f\{E) involved in the decompo- 



sition (|131|) reads in terms of Bessel functions Jy and Yy 
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Yiikuo) J2[ku{E)] - Ji{kuo) Y2[kuiE)] 



(136) 



with k 



A-l 
D 



and 



A4 



and (|131D may be written with the notation = 1 + Dk^ as 

POO 

P{E, t) = p{E) + 2De-^°* / dk a(Afc) /a, (^) e 



(137) 



;i38) 



In the hmit — > 0, the behaviors /a^. oc k"^ and a(Afc) oc k'^ for a generic initial condition (see 
120| ) gives the asymptotic algebraic correction to the exponential e~^°^ 



P{E,t)-p{E)^ -—- 



(139) 



to be compared with (|^). The presence of the diffusion term in (0) representing the interaction 
between particles thus tends to accelerate the relaxation in comparison to the case (|I]) of 
independent particles. 
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Figure Captions 

Fig 1 Correlation function Ccq(t) (multiplied by a 'mock' /3 relaxation C/3(t) = exp — ^-^j^, 
where r{t) describes a diffusive motion in an harmonic potential well: r^(t) = ^q[1 — exp(— ;^)]; 
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^0 can be thought as the 'size' of the cage, and ^ of the order of the high temperature diffusion 
constant) . 

Fig 1-a Exponential density of states. 

Plot of CeqC^(t) versus log]^o(ro^) for q^Q = 0.5, Tq = SFq ^, and ^ = 2., 1.1, 1.03. 
Fig 1-b Gaussian density of states. 

Plot of CeqC/3(t) versus logio(rot) for q^o = 0.5, tq = 5ro\ and |^ = 0.5, 0.25, 0.15. Note 
that the plateau observed for the lowest temperature eventually decays to zero. 
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